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Preface 

The  altimeter  described  within  necessarily  uses  a certain  sensor  to  obtain  the 
results  reported.  However,  this  report  is  not  endorsing  that  sensor  but  rather 
describing  a means  of  using  any  sensor  with  a similar  output  signal. 
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contributions  noted.  He  also  is  indebted  to  Mr.  R.  Cowie,  Lt.  L.  Wrinkle,  and 
Mrs.  C.  Rice  for  their  constructive  remarks  and  suggestions  on  the  first  draft; 
and  to  Ms.  M.  Cross  for  typing  the  manuscript. 
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An  Autoranging  Balloon  Altimeter: 
A Single  Pressure  Transducer 
Monitors  Altitude  From  0 to  44  Kilometers 
With  30  Meters  Resolution 


I.  IM  HOm  CTION 


The  purpose  of  this  report  is  to  document  the  use  of  a single  pressure  sensor 
as  an  altimeter  for  balloon  borne  scientific  experiments  in  the  altitude  range  0 to 
44  m (145  k ft).  The  report  describes  the  development  of  an  algorithm  to  facilitate 
use  of  the  sensor,  and  the  fabrication  of  an  interface  device  to  perform  the  algo- 
rithm. The  sensor  and  interface  operate  into  an  existing  frequency-sensitive  en- 
coder with  a commutator  segment  aperture  time  of  2 sec. 


2.  Till:  SKNSOII 


The  Model-ID  Digital  Pressure  Transducer  by  Hamilton  Standard^  is  a small, 
lightweight  sensor  with  a variable  frequency  output.  The  calibration  error  is 
±0.  008  percent  of  full  scale  (FS)  and  the  repeatability  error  is  ±0.  0001  percent  FS. 
Since  full  scale  pressure  is  20.0  psia,  the  absolute  calibration  error,  E,  is  0.0032 
psia.  To  relate  this  pressure  error  to  altitude  error,  we  chose  to  define  the  de- 
pendence of  atmospheric  pressure  upon  altitude  by 


(Received  for  publication  24  January  1978) 

1.  Hamilton  Standard  Instruction  Manual:  The  Hamilton  Standard  Model  ID  Digital 
Pressure  Transducer  (P/N  752200),  Windsor  Locks,  Connecticut. 
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HVSEDINa  PifSl  BUlNK-WOT  FI, 


(1) 


P = P e°^ 
o 

where  P is  pressure  in  psia,  P^  is  14.  696  psia,  cc  is  a constant  equal  to 

-0.  04401 1/k  ft  when  A is  altitude  in  k ft.  This  is  a very  good  approximation  of  the 

2 

atmosphere  as  defined  in  the  U.S.  Standard  Atmosphere  1976.  Eq.  (1)  is  much 
easier  to  work  with  than  the  altitude  "shells"  defined  in  reference  2,  computational 
equations  section.  Differentiating  Eq.  (1)  yields 


r 


dP 

dA 


(2) 


where  r is  the  rate  of  change  in  psia/ft  and  r equals  -0.  000646786  psia/ft.  There- 
fore, the  absolute  error  in  feet,  E',  is  defined  as 

E'  = E/r.  (3) 

Table  1 was  generated  using  Eq.  (3),  and  shows  that  this  sensor  fits  our  re- 
quirements very  well.  However,  Table  2 and  Figure  1 show  that  the  output  rate 
of  change  decreases  rapidly  with  increasing  altitude.  Therefore,  some  manipula- 
tion is  required  to  maximize  the  sensor's  utility. 


Table  1.  Sensor  Error 


Altitude 

Absolute  Error 

km 

(k  ft) 

ft 

m 

0 

0 

4.  95 

1.  51 

6.  1 

20 

11.  9 

3.  64 

12.  2 

40 

28.  8 

8.  77 

18.3 

60 

69.4 

21.  1 

24.4 

80 

167 

51.  0 

30.  5 

100 

403 

123 

36.  6 

120 

973 

297 

42.  7 

140 

2346 

715 

2.  U.S.  Standard  Atmosphere  (1976).  NOAA,  NASA,  USAF,  Washington,  D.C. 
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Figure  1.  Sensor  Output:  Period  vs  Altitude 


The  manufacturer  defines  the  sensor  in  terms  of  two  equations.  The  first 
equation 

Y = B + B,  X , 
o 1 

produces  a dummy  variable  Y as  a function  of  x,  the  output  period  in  microseconds. 
The  second  equation 

P = A^  + Aj  Y + A2  Y^  + A3  Y^  + A^  Y^ 

defines  the  pressure  (psia)  as  a function  of  the  dummy  variable.  Table  2 was  gen- 
erated using  the  vender's  tables  for  the  sensor  at  hand  and  pressures  from  the 
U.  S.  Standard  Atmosphere  1976.  Periods  were  calculated  via  linear  interpolation 
for  those  pressure  values. 
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Table  2.  Sensor  Output  (Preliminary) 


Altitude 

Pressure 

psia 

Period 

^ sec 

km 

k ft 

-0.  30 

-1 

15.  234 

194. 954 

0.  00 

0 

14. 696 

196.  043 

1.  52 

5 

12.  228 

201.329 

3.  05 

10 

10.  108 

206.  298 

4.  57 

15 

8.  297 

210.911 

6.  10 

20 

6.  758 

215.  140 

7.  62 

25 

5.  460 

218. 964 

9.  14 

30 

4.  372 

222.374 

10.  7 

35 

3.467 

225.373 

12.  2 

40 

2.  730 

227. 939 

13.  7 

45 

2.  148 

230.052 

15.  2 

50 

1.  691 

231. 770 

16.  8 

55 

1. 331 

233. 165 

18.  3 

60 

1.  048 

234. 288 

19.  8 

65 

0.  825 

235.  188 

21.3 

70 

0.  650 

235.  908 

22.  9 

75 

0.  512 

236. 479 

24.  4 

80 

0.  403 

236.  938 

25.  9 

85 

0.  318 

237.300 

27.  4 

90 

0.  251 

237.  583 

29.  0 

95 

0.  200 

237.  804 

30.  5 

100 

0.  160 

237.974 

32.  0 

105 

0.  128 

238.  113 

33.  5 

1 10 

0.  103 

238. 222 

35.  1 

115 

0.  0840 

238. 304 

3 6.  6 

120 

0.  0683 

238.373 

38.  1 

125 

0.  0557 

238. 427 

39.  6 

130 

0.  0456 

238. 471 

41.  1 

135 

0.  0375 

238. 507 

42.  7 

140 

0.  0309 

238. 535 

44.  2 

145 

0.  0255 

238. 559 
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3.  THE  ALGORITHM 


Prior  to  discussing  the  proposed  algorithm,  two  constraints  must  be  consid- 
ered. First,  the  altitude  resolution  corresponding  to  ±1/2  bit  error  of  the  binary 
interface,  should  be  approximately  30  m (100  ft)  worse  case.  Second,  the  data 
must  fit  the  twelve  bit  register  in  the  aforementioned  encoder.  These  twelve  bits 
correspond  to  4096  states  or  eight  512  state  ranges.  The  reason  for  this  will  be- 
come evident  as  the  discussion  continues. 

It  was  decided  to  multiply  and  measure  the  period  rather  than  multiply  and 
measure  the  frequency.  This  approach  was  taken  because  it  has  two  strong  posi- 
tive points.  First,  it  is  easier  to  accomplish  with  the  technology  chosen  for  device 
fabrication.  Second,  it  is  easier  to  achieve  a uniform  pulse  rate  which  facilitates 
a uniform  predictable  error  (resolution)  at  any  given  sensor  frequency. 

In  general  terms,  the  algorithm  is  based  on  the  fact  that  an  unknown  period 
can  be  measured  by  gating  a known  reference  frequency  by  that  period  while  count- 
ing the  number  of  events  which  occur.  The  resulting  count  is  a measure  of  the 
period  with  an  error  of  one  event.  Therefore,  enlarging  (by  multiplying)  the  un- 
known period  to  increase  the  number  of  events  of  the  reference  oscillator  occurring 
in  the  unknown  time,  produces  a large  count  which  reduces  the  relative  measure- 
ment error.  The  only  drawback  with  this  scheme  is  that  the  large  constant  portion 
of  the  sensor  period  (approximately  190  ^sec)  is  multiplied  along  with  the  changing 
period  (approximately  50  /xsec)  which  contains  the  information.  Rather  than  trans- 
mit this  unchanging,  known  and  therefore  informationless  number  of  events,  it  will 
be  removed  by  subtraction.  Figure  2 is  a block  diagram  of  the  algorithm. 

The  principal  functions  are  in  line  under  the  sensor  block;  multiply  the  output 
period  (by  2*^),  measure  the  resulting  period,  and  subtract  the  known  part  of  the 
count.  After  subtraction,  the  remainder  will  be  less  than  512  events  of  altitude 
information  for  any  scale.  As  indicated  earlier,  4096  events  are  eight  512  event 
subsets;  and  as  may  be  surmised  from  Figure  2,  one  coefficient  for  the  period 
multiplier  will  not  suffice,  nor  will  one  subtrahend  for  the  subtractor.  Therefore, 
a mechanism  to  change  these  constants  (that  is,  constants  within  each  range)  must 
be  supplied;  this  is  the  scale  selector. 

Without  considering  how  it  will  be  accomplished.  Table  3 lists  events  vs  alti- 
tude for  each  scale.  The  multiplier  2^  X 10®  is  a combination  of  the  period  multi- 
plier 2*^  (unitless)  with  the  product  modulator  10®  (events/sec). 

Note  that  the  resolution  varies  from  about  17  m (55  ft)  to  32  m (105  ft)  per 
scale;  and  there  is  plenty  of  overlap  in  the  ranges.  By  this,  I mean  that  the  "6" 
scale  could  be  used  for  altitudes  over  130  k ft  even  though  the  instrument  is  sup- 
posed to  be  switched  to  the  ”7"  scale  by  that  point.  And,  the  "7"  scale  is  func- 
tional below  120  k ft  event  through  it  is  supposed  to  be  used  from  125  to  145  k ft. 
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Table  3.  Data  Summary  (Preliminary) 
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Table  3.  Data  Summary  (Preliminary)  (Cent.) 


Scale 

alt 

p X multiplier 

= events 

A 

Resolution 

Remaining 

Counts 

k ft 

^ sec  2"  X lOfi 

events 

ft 

M 

243  939 

59 

139 

6 

BB 

n = 10 

244  023 

70 

223 

subtrahend:  243, 800 

244  093 

55 

90 

293 

IM 

244  149 

45 

110 

349 

130 

244  194 

394 

120 

488  187 

111 

45 

87 

125 

n = 11 

488  298 

90 

55 

198 

130 

subtrahend:  488,  100 

488  388 

73 

67 

288 

■H 

135 

488  462 

57 

87 

362 

140 

488  519 

49 

101 

419 

145 

488  586 

468 

After  calculation  of  the  events  and  the  resolution,  the  subtrahend  is  chosen  to 
place  the  counts  remaining  in  the  range  of  0 to  512  for  each  scale. 

How  does  the  scale  selector  work  and  what  is  its  switching  resolution?  Fig- 
ure 2 shows  that  the  scale  selector  has  inputs  from  the  period  multiplier  and  the 
reference  oscillator.  The  product  (sensor  period  x multiplier  x coefficient)  is 
formed  and  then  compared  to  several  predetermined  products  which  are  defined  as 
switching  points.  The  most  significant  point  to  be  detected  determines  the  scale 
to  be  used.  Table  4 lists  the  switching  points  and  resolution.  The  resolution  was 
computed  by  determining  the  data  count  for  a point  on  either  side  of  the  switching 
point,  finding  the  difference,  and  dividing  the  altitude  difference  by  the  count  dif- 
ference. A worse  case  example  shows  that  in  the  vicinity  of  a switching  point, 
usable  data  are  available  on  both  scales.  Note  that  at  125  k ft,  the  resolution  is 
304.8  m (1000  ft)*:  on  Table  3 it  can  be  seen  that  this  corresponds  to  9 counts 
above  125  k ft  in  scale  "6"  or  22  counts  below  125  k ft  in  scale  "7".  There  is 
plenty  of  room  in  the  event  counter  to  accommodate  this  resolution. 

One  loose  end,  differentiating  between  scales,  remains.  Recall  that  the  data 
have  been  confined  to  512  events,  while  the  data  counter  in  the  encoder  accommo- 
dates 4096  (8  X 512)  states.  The  data  counter's  output  are  a twelve  binary  bit  word 
which  is  arranged  in  four,  3 bit  groups  to  facilitate  recording  and  decoding.  Nine 
bits  define  512  states  (0  through  511  events).  Since  the  data  remain  below  512 

This  is  an  exact  definition  from  reference  2. 
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Table  4.  Scale  Selector  Switching  Point  Resolution 


Altitude 

Pressure 

Period 

Multiplier 

Resolution 

k ft 

psla 

4 sec 

A 

ft 

30 

29.  5 

4.  472 

222.  052 

28422 

82 

12 

30.  5 

4.  274 

222.  690 

28504 

50 

49.  5 

1. 732 

231. 615 

29646 

40 

25 

50.  5 

1.  651 

231. 923 

29686 

65 

64.  5 

0.  845 

235.  108 

30093 

13 

48 

65.  5 

0.  806 

235.  267 

30114 

80 

79.  5 

0.413 

236.  900 

30323 

9 

111 

80.  5 

0.394 

236.  975 

30332 

95 

94.  5 

0.  205 

237. 781 

30435 

5 

200 

95.  5 

0.  196 

237.  819 

30440 

no 

109.  0 

0.  108 

238.  200 

30489 

5 

400 

111.0 

0.  0995 

238.  237 

30494 

125 

124.  0 

0.  0580 

238.  417 

30517 

2 

1000 

126.  0 

0.  0535 

238. 437 

30519 

events,  the  most  significant  three  bits  are  never  used:  or,  they  always  read  000. 

If  512  events  were  added  to  the  data,  the  most  significant  three  bits  would  contain 
001.  Note  that  this  set  of  bits  has  now  uniquely  defined  two  scales:  a 0 scale,  and 
a 1 scale.  If  integer  multiples  of  512  (1024,  1536,  2048,  etc.)  are  added  to  the 
data,  more  scales  (2,  3,  4,  etc.)  are  uniquely  defined.  This  method  defines  eight 
scales  designated  000  to  111  in  the  binary  data  register  under  discussion. 

Defining  new  subtrahends  implements  this  method  of  scale  definition  in  the 
sensor  interface  algorithm.  Since,  events  output  = total  count  - subtrahend,  we 
add  512  counts  to  the  output  for  scale  1.  Then  events  output  + 512  = total  count  - 
subtrahend  + 512  or  events  output  ♦ 512  = total  count  - [subtrahend  - 512)  = total 
count  - new  subtrahend  where  the  "new  subtrahend"  is  512  counts  less  than  the 
original  subtrahend.  Likewise,  by  removing  multiples  of  512  from  the  original 
subtrahends  defined  in  Table  2,  new  subtrahends  are  defined  which  also  act  as 
scale  identifiers  (see  Table  5).  The  resolution  per  scale  remains  as  listed  in 
Table  3. 
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Table  5.  Output  Data  Summary 


Scale 

Altitude 

Preliminary 

Subtrahend 

New 

Subtrahend 

Counts 

unitless 

km 

k ft 

unitless 

0 

0.30 

-1 

3100 

3100 

19 

0.00 

0 

36 

1.  52 

5 

121 

7.  62 

25 

401 

9.  14 

30 

457 

1 

9.  14 

30 

7000 

6488 

627 

10.  7 

35 

723 

12.  2 

40 

806 

13.7 

45 

873 

15.  2 

50 

928 

2 

15.  2 

50 

14750 

13726 

1107 

16.  8 

55 

1196 

18.3 

60 

1268 

19.  8 

65 

1326 

3 

19.  8 

65 

30000 

28464 

1640 

21.3 

70 

1732 

22.  9 

75 

1805 

24.  4 

80 

1864 

4 

24.  4 

80 

60500 

58452 

2204 

25.  9 

85 

2296 

27.  4 

90 

2369 

29.0 

95 

2425 

5 

29.0 

95 

121600 

119040 

2715 

30.  5 

100 

2802 

32.  0 

105 

2873 

33.  5 

110 

2929 

r> 

33.  5 

110 

243800 

240728 

3211 

35.  1 

115 

3295 

36.  6 

120 

3365 

38.  1 

125 

3421 

7 

38.  1 

125 

488100 

484516 

3782 

39.  6 

130 

3872 

41.  1 

135 

3946 

42.  7 

140 

4003 

44.  2 

145 

4052 

IB 


Note  that  Tables  3 and  5 were  begun  with  an  altitude  below  sea  level  to  in- 
sure that  the  normal  daily  variation  in  barometric  pressure  at  sea  level  locations, 
would  not  underflow  the  data  register  and  produce  confusing  results.  The  output 
data  summarized  in  Table  5 is  graphed  in  Figure  3. 


Figure  3.  Count  vs  Altitude 


1.  H X TIONAI,  ANALYSIS 
1. 1 lieneral 

With  the  description  of  the  algorithm  behind  us,  it  is  time  to  consider  each 
of  the  functional  blocks  of  Figure  2.  At  this  point  they  will  not  be  analyzed  to  the 
point  of  what  device  implements  what  function,  but  rather  to  the  level  of  logic  func- 
tions needed  to  accomplish  each  task.  The  device  number  and  pin  utilization  will 
be  considered  in  a later  section.  Positive  logic  will  be  used  throughout. 

Figure  4 is  based  on  Figure  2 and  considers  power  supplies,  and  the  number 
of  wires  to  get  information  from  block  to  block.  Power  supply  wires  and  the 
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SYNC 

OUT 


Figure  4.  Equipment  Block  Diagram 


initializer  pulse  routing  are  not  shown;  consider  that  when  they  are  needed,  they 
will  be  available.  The  capital  letters  are  placed  at  break  points  for  referencing. 

1.2  \ <illa;'c  lli‘;;uliitors 

All  power  will  be  supplied  from  a 12  Vdc  battery  with  negative  system  ground. 
Sensor  power,  ±15  Vdc,  will  be  obtained  from  a commercially  available  dc/dc  con- 
verter. I.ogic  power,  1 10  Vdc,  will  be  supplied  by  a regulator  documented  in  the 
next  section. 

1.2  Initiidi/.iT 

A logic  "1"  will  be  generated  by  an  RC  differentiator  when  power  is  applied 
and  will  be  used  where  necessary  to  insure  that  the  state  of  certain  devices  is 
known  when  the  measuring  process  begins. 

t.  i HffiTcncr  (I'M'!  Iliilor 

A one  megahertz,  crystal  controlled,  capacitor  trimmable  oscillator  with 
50-50  duty  cycle  for  CMClS  logic  is  commercially  available.  Utilized  in  a voltage 
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regulated,  constant  load  regime  with  small  temperature  change  (20°C  ± 10°C),  a 
maximum  error  of  ±1/2  PPM  is  possible. 


1.5  Sensor 

■Several  characteristics  of  the  sensor  have  been  covered  in  the  preceding  sec- 
tions and  more  information  has  been  referenced.  1 he  most  significant  point  now 
is  that  the  output  swings  from  0 to  5 Vdc  in  a 50-50  duty  cycle  square  wave.  This 
voltage  level  must  be  up-shifted  to  the  10  Vdc  level  at  which  the  CMOS  logic 
operates. 


1.6  Level  -Shifler 

A simple  transistor  inverter  will  perform  the  level  shifting  operation.  It  will 
be  followed  by  an  inverting  buffer  to  insure  sharp  transitions  (see  Figure  5), 

L7  Period  Multiplier 

This  function  is  performed  by  a synchronous  binary  counter  (see  Figure  6). 
Note  that  twelve  stages  of  period  multiplication  are  shown,  even  though  Table  3 
indicates  that  eleven  is  the  highest  power  used.  The  reason  for  this  is  shown  in 
Figure  6(b).  Carrying  the  multiplication  one  step  further  than  necessary  results 
in  a series  of  logic  " 1”  (or  "0")  pulses  which  has  the  duration  of  the  period  of  the 
required  state.  The  utility  of  this  approach  will  be  shown  in  a few  paragraphs. 
Note  that  the  multiplier  is  initialized  when  power  is  applied  and  then  tracks  the 
period  of  the  sensor.  The  output  of  any  Q is 

Qn  PX2"  , 

where  p is  the  period  of  the  sensor;  the  units  of  Q are  the  same  as  p. 

I.K  Prodiicl  Modulator 

This  section  includes  a commutator  and  multiplier.  As  seen  in  Figure  7,  the 
multiplier  is  a simple  AND  gate  which  by  definition  performs  the  multiplication 
function  as  indicated  by  its  symbols;  dot  (. ) and  cross  (x).  Another  way  of  looking 
at  the  process  stems  from  an  older  name  of  the  AND  circuit;  a coincidence  cir- 
cult.  When  both  inputs  are  logic  I's,  a logic  1 appears  at  the  output;  or,  the 
input  modulates  the  reference  oscillator.  Hence,  the  name  describes  the  func- 
tion; (iroduct  modulator. 

•See  reference  3,  page  317. 

3.  Millman,  J.,  and  Taub,  H.  (1965)  Pulse,  Digital,  and  .Switching  Waveforms, 
McGraw-Hill  Book  Co.,  New  York. 
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Tri-state  transmission  gates  form  the  commutator  which  is  controlled  by  the 
scale  selector.  The  output  E is  used  as  a flag  to  insure  that  the  scale  selector 
does  not  change  commutator  segments  while  a multiplication  is  being  performed. 

1.')  Scale  Selector 

As  indicated  in  Figure  2,  the  scale  selector  monitors  the  period  multiplier 
and  the  reference  oscillator:  it  then  performs  a period  measurement  through  the 
product  multiplier.  Since  the  period  multiplier  output  (Qg)  and  reference  oscilla- 
tor output  are  constants,  the  resolution  of  the  selector  tracks  the  frequency  out- 
put from  the  sensor  and  decreases  with  increasing  altitudes.  This  was  defined  in 
Table  4.  A read  only  memory,  TOM,  is  used  to  detect  the  proper  set  of  outputs 
which  defines  the  number  of  events  in  Table  3.  There  is  one  ROM  output  for  each 
of  the  seven  switching  points.  An  R/S  latch  is  used  to  detect  a logic  1 at  the  ROM 
outputs  (see  Figure  8).  A bank  of  gates  allows  only  the  most  significant  (highest 
magnification)  scale  detected  to  be  presented  to  the  clocked  latches  which  form  the 
output  of  the  scale  selector.  After  a scale  has  been  selected  and  clocked  into  the 
D latches,  the  counter  and  R/S  latches  are  reset  for  the  next  selection  process. 

Figure  8 shows  that  six  outputs  are  derived  from  the  Q outputs  brought  in  to 
the  scale  selector  as  group  C:  Figure  9 shows  where  they  fall  in  relation  to  the 
period  being  measured  (shaded  Qg  pulse).  Note  that  the  working  (VVRK)  period  is 
not  used  directly,  but  defines  the  time  during  which  the  jam  (JAM)  and  reset  (RST) 
pulses  are  generated.  Remerriber  that  the  Q outputs  are  wholly  dependent  on  the 
period  of  the  sensor,  thereby  locking  the  JAM  and  RST  pulses  to  the  same  stand- 
ard. Therefore,  the  scale  is  being  selected  continuously,  based  on  the  frequency 
of  the  sensor;  meanwhile,  the  sensor  period  is  being  measured  for  the  output 
(point  G)  by  another  part  of  the  instrument.  The  scale  cannot  be  changed  during  a 
measurement  cycle,  defined  by  E (see  Figures  G and  7),  because  the  signal  at  E 
is  used  to  inhibit  the  JAM  pulse  when  the  product  multiplier  and  subtractor  are 
busy. 

t.lO  Sut)lrd<  li)r* 

The  subtracter  and  multiplier  work  together  to  produce  a pulse  burst  as  pre- 
viously described.  Figure  10  has  an  expanded  block  diagram  of  the  subtracter 
followed  by  a simple  block  schematic.  Besides  subtracting  the  requisite  number 
of  events  from  the  pulse  burst  produced  by  the  product  modulator,  the  subtracter 
also  synchronizes  the  output  with  the  encoder.  Therefore,  the  latter  part  of  this 
subsection  is  heavily  slanted  toward  presenting  the  information  to  that  encoder. 


See  Appendix  D. 
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a)  EXPANDED  BLOCK  DIAGRAM 


03  iumr ■ ••■juuuL  - 


WRK  I I I L 
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n 

Figxire  9.  Scale  Selector  Timing 


In  another  scheme  the  informative  number  of  events  could  very  easily  be  stored  in 
the  counters  for  subsequent  processing.  Figure  11  presents  waveforms  that  aid 
understanding  the  subtracter's  functioning. 

As  explained  earlier,  when  a certain  output  is  of  interest,  the  output 
is  used  to  facilitate  mixing  by  the  AND  circuit.  Figure  11  has  been  scaled  for  a 
period  of  about  half  a second,  which  is  slightly  larger  than  that  produced  in 
scale  "7"  operation.  The  J input  (logic  1)  is  2 sec  long.  The  three  waveforms  in 
part  (a)  are,  as  explained  in  previous  section,  keyed  to  the  sensor  and  continuing. 
The  first  waveform  is  the  product  modulator  output:  the  minuend  M where 
M Qfj+l  ’ ^ coefficient  C is  10*’  events/second. 

The  subtractor  is  controlled  by  the  J input  from  the  encoder  through  the  per- 
iod-detecting dual  flip-flop  circuit.  When  gate  1 is  enabled,  the  counter  presents 
to  the  ROM  the  elapsed  events.  The  scale  selector  chooses  the  proper  ROM  sub- 
trahend. When  that  number  appears,  (waveform  h)  the  latch  opens  gate  2,  pro- 
ducing waveform  G which  is  a pulse  burst  containing  the  difference  between  the 
minuend  and  the  subtrahend.  Figure  11(b)  and  (c)  show  the  effect  of  changing  the 
timing  of  J.  In  each  case,  the  difference  is  the  same,  but  the  placement  of  the 
pulse  burst  within  the  encoder  aperture  (j)  changes. 


a)  EXPANDED  BLOCK  DIAGRAM 


b)  LOGIC  DIAGRAM 


Figure  10.  .Subtracter  (Basic) 
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5.  IMPI.EMKNTMION 


5.1  (JeinTal 

This  section  explains  a collection  of  schematics  suitable  for  fabrication;  that 
is,  it  gets  down  to  the  package  level.  The  schematics  are  in  Appendix  A.  Most 
of  the  circuits  are  straight  forward;  and,  with  the  explanation  of  section  4,  should 
be  understandable.  Therefore,  unless  a point  needs  further  clarification,  it  will 
not  be  discussed. 

5.2  Initialization 

The  schematic  block  diagram  shows  that  the  reset  pulse  generated  by  the 
initializer  goes  only  to  the  period  multiplier.  Why?  How  do  the  counters  in  the 
scale  selector  and  subtractor  get  initialized?  The  reset  pulse  (RST)  generated  by 
the  scale  selector  control  section  initializes  them. 

The  initializer  holds  the  period  multiplier  reset  for  about  12  msec.  Subse- 
quently, the  counter  begins  to  advance;  the  WRK  output  is  a logic  1 because  suffi- 
cient events  have  not  occurred.  In  2.  4 msec,  an  0.  8 msec  RST  pulse  initializes 
the  scale  selector  and  subtractor.  These  times  are  based  on  a 0.  20  msec  sensor 
period.  RST  pulses  occur  every  51  msec,  so  the  scale  data  will  be  fresh  by  time 
the  encoder  completes  its  initialization  and  framing  pulse  generation  (2  sec).  See 
Figure  12.  Note  that  these  times  are  based  on  a specific  sensor  output  period 
which  was  assumed  to  have  a logic  1 to  0 transition  immediately  before  the  initial- 
izing pulse  ended, 

5.5  ('.apacitors 

Every  capacitor  in  this  device  is  for  power  supply  filtering  or  radio  frequency 
energy  (RF)  bypassing.  The  latter  are  included  because  of  the  severe  RF  environ- 
ment in  which  the  altimeter  must  function. 

5.  t l.opir  .Sub.'.ilitutions 

In  section  4,  it  was  convenient  to  define  logic  functions  with  AND  gates  and 
NOT  gates.  However,  during  implementation  this  is  usually  not  optimum.  For 
example,  the  output  at  point  F in  Figure  A4  is  inverted  with  respect  to  an  AND 
output;  but,  since  there  is  no  attempt  to  control  the  phase  relationship  of  Y1  and 
the  sensor,  it  doesn't  matter.  In  other  cases,  the  need  for  an  AND  gate  was  im- 
perative so  it  was  provided  as  in  Figure  A 5.  In  this  same  figure,  note  that  the 
NOR  gate  is  the  implementation  of  the  AND  gate  with  inverted  inputs  in  Figure  8b. 
Such  conversions  are  in  accordance  with  DeMorgan's  laws. 

*See  reference  3,  page  328-330. 
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Figure  12.  Interface  Start-up  Waveforms 

ROM  I’ro^rainmin)' 

There  are  two  ROM's  in  this  instrument:  the  scale  selector  ROM  will  be  dis- 
cussed first.  This  ROM,  which  detects  when  the  count  for  each  of  the  scale  "break 
points"  is  reached,  is  programmed  for  the  binary  count,  as  noted  in  Table  (;(a). 

The  ROM  passes  the  information  from  the  counter  to  the  latches  virtually  instan- 
taneously; the  latch  reacts  in  nsec. 

This  highlights  a potential  problem  in  the  Subtractor  ROM  function.  If  that 
ROM  were  set  to  detect  the  subtrahends  listed  in  Table  4,  the  event  at  which  the 
subtrahend  was  detected  would  be  counted  as  part  of  the  remainder.  This  is  pre- 
cluded by  setting  the  ROM  for  the  number  equal  to  the  subtrahend  plus  1.  These 
values  and  the  resulting  codes  are  listed  in  Table  7(a).  The  second  part  of  each 
Table  (6  and  7)  lists  the  pin  number  links  to  be  removed  from  each  new  matrix  to 
produce  the  defined  ROM.  These  same  links  are  redefined  by  row  and  column  in 
Table  8. 


30 


Table  6.  Scale  Selector  ROM  Codes 


(a)  binary  word  generation 


Break  Points 

Period 

sec 

Multiplier 

Hertz 

Count 

k ft 

km 

events 

eventSg 

events  2 

30 

9.  14 

222.  374 

28463 

67457 

110  111  100  101  111 

50 

15.  2 

231. 770 

|_ 

29666 

71742 

111  001  111  100  010 

65 

19.  8 

235.  188 

30104 

72630 

111  010  1 10  on  000 

80 

24.4 

236.  938 

30328 

73170 

111  on  001  111  000 

95 

29.0 

237. 804 

30438 

73346 

111  on  on  100  no 

110 

33.  5 

238. 222 

30492 

73434 

111  on  100  on  loo 

125 

38.  1 

238. 427 

mM 

30518 

73466 

111  on  100  no  no 

(b)  paths  to  remove  (pin  numbers) 


Column 

A12-A14 

Rows 

A12 

A13 

A14 

1 

6 

3,4 

4 

7 

2,  4,  5,  6 

none 

2,  3 

14 

2,  3,  4 

2,  3.  6 

3 

8 

2,  3,  4 

4,  5 

3 

13 

2.  5,  6 

5 

3 

9 

2,  3 

2.3.4 

3 

12 

2,  5 

3,4 

3 
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Table  7.  Subtracter  ROM  Codes 


(a)  binary  word  generation 


HWi 

Subtrahend  +1 

mM 

base  10 

base  8 

base  2 

0 

3101 

6035 

000  000  000  110  000  on  loi 

1 

6489 

14531 

000  000  001  100  101  on  ooi 

2 

13727 

32637 

000  000  on  010  noon  in 

3 

28465 

67461 

000  000  no  111  100  no  ooi 

4 

58453 

162125 

000  001  no  010  001  oio  loi 

5 

119041 

350401 

000  on  101  000  100  ooo  ooi 

6 

240729 

726131 

000  ill  010  no  001  on  ooi 

7 

484517 

1662245 

001  no  no  oio  oio  loo  loi 

(b)  paths  to  remove  (pin  numbers) 


Column 

A26-A29 

Rows 

A26 

A27 

A28 

A 29 

1 

3 

2,  3,  4,  5,  6 

4,  5,  6 

2,3,4.  5 

7 

3,  4 

2,  4,  6 

2,  5,  6 

2,  3,  4,  5 

14 

none 

2,  3,  6 

3.  6 

2,  3,  4,  5 

8 

3,  4,  5 

3.4 

4 

2,  3.  4,  5 

13 

3.5 

2,  4,  5,  6 

3,  4 

3,  4,  5 

9 

3,  4,  5,  6 

2,  3,  4,  6 

2,  3,  5 

4,5 

12 

3,4 

2,  4,  5,  6 

4,6 

5 

10 

3,  5,  6 

3,  5,  6 

3,  4 

2 

From  Table  4. 
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Table  8.  Scale  Selector  and  Subtractor  ROM  Paths  Removed  by  Row  and  Column 


(a) 


(b) 


Column 

A12-A14 

Rows 

A12 

A13 

A14 

1 

5 

2,3 

3 

2 

1,3,4,  5 

^one 

1,2 

3 

1,2,3 

1,2,  5 

2 

4 

1,2,3 

3,4 

2 

5 

1,  4,  5 

4 

2 

6 

1,2 

1,2,3 

2 

7 

1,4 

2,  3 

2 

.l.h  h aliriratiun 

The  transducer  and  the  circuits  documented  in  Appendix  A are  packaged  in  a 
If).  5 cm  X 14.  0 cm  X 10.  2 cm  (G.  5 in  X 5.  5 in  X 4.  0 in)  container.  The  circuts  are 
on  three  circuit  cards  whose  layouts  are  in  Appendix  B;  the  instrument  is  shown  in 
Figure  13. 
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As  described  earlier,  the  data  word  has  twelve  bits  in  four  3 -bit  groups.  The 
existing  data  acquisition  system,  into  which  the  sensor  and  its  interface  gain  access 
through  the  encoder,  prints  data  words  in  a letter  code  defined  by  Table  9;  each 
data  word  is  printed  as  a four  letter  code.  The  object  of  the  data  retrieval  system 
is  to  relate  every  possible  code  word  to  the  U.  S.  Standard  Atmosphere;  a computer 
program  accomplished  this  task. 


Table  9.  Codes 


Binary 

Octal 

Letter 

000 

0 

S 

001 

1 

U 

010 

2 

R 

Oil 

3 

W 

100 

4 

D 

101 

5 

K 

no 

6 

G 



7 

O 

Figure  14  is  the  program  flow  chart.  Part  (a)  is  the  overall  program  while 
part  (b)  details  the  atmospheric  shell  determination.  Table  10  contains  the  ref- 
erence  constants  used  to  define  the  shells  and  computes  the  exact  altitude.  Table 
11  defines  the  mnemonic  codes  used  in  the  program. 

Compare  the  flowchart  to  the  sensor  interface  algorithm  and  note  that  it  is  the 
inverse  of  the  algorithm.  The  outside  (largest)  DO  loop  defines  the  eight  scales: 
the  next  four  DO  loops  define  the  512  bits  in  each  scale.  In  this  way  every  possible 
code  combination  is  listed  even  though  many  can  only  be  generated  if  the  scale 
selector  is  disconnected,  overridden,  or  malfunctioning.  The  51 2 bits/scale  are 
generated  by  four  DO  loops  to  facilitate  a neat  printout.  The  source  code  and  sub- 
sequent machine  listing  are  contained  in  Appendix  C. 

The  body  of  the  lookup  table  or  dictionary  is  less  than  8 in.  X 10  in.  so  it  can 
be  cut  to  fit  a standard  ring  binder.  Average  line  resolution  data  is  on  the  far 
right. 


Based  on  U.  .S. 


Standard  Atmosphere, 


reference  2. 
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Table  10.  Atmospheric  Constants 


Subscript 

Geopotential 

Height 

T emperature 
Gradient 

Pressure 

T emperature 

b 

»b 

^m.  b 

Pb 

T u 
m.  b 

unitless 

km' 

k ft 

K/km' 

mb 

K 

0 

0 

0 

-6.  5 

1013. 2500 

288. 15000 

1 

11 

36. 089239 

0.  0* 

226. 321 

216.  6500 

2 

20 

65.  616798 

+ 1.  0 

54. 7489 

216.  6500 

3 

32 

104. 98688 

+2.  8 

8. 68014 

228. 6500 

4 

47 

154. 19948 

0.  0 

1.  1090 

270. 6500 

4 

In  the  computer  program,  0.  0001  in  place  of  0.  0 yields  data  identical  to  that 
tabulated  in  reference  2 and  precludes  using  a different  equation.  The  author 
thanks  Mr.  James  F.  Dwyer  for  this  suggestion. 


Table  11.  Mnemonic  Codes 


Code 

Definition 

IBK 

block,  8 lines  of  code 

IL 

leaf  or  page,  four  per  scale 

IPG 

page  counter 

I SB 

subtrahend 

LSC 

scale,  MSB,  eight  per  sensor 

ITC 

total  count,  same  as  minuend 

J2 

block  counter,  2 per  page 

J3 

line  counter,  8 per  block 

J4 

column  counter,  LSB,  8 per  line 

PR 

pressure,  atmospheric 
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7.1  Kum  tional  Tests 

Prior  to  checking  how  well  (or  poorly)  the  device  performed,  it  was  necessary 
to  see  if  it  worked  at  all.  Table  12  is  the  general  procedure  used  to  check  the 
instrument.  The  first  step  is  critical;  there  is  no  margin  for  error. 


Table  12.  Check  Procedure 


.Step 

Function 

1 

With  power  disconnected  assure  that  the  sensor 
mating  connector  is  disconnected.  Apply  power, 
measure  voltage,  and  check  polarity  at  the 
sensor  mating  connector. 

2 

Check  initializer  pulse  at  period  multiplier. 

3 

Connect  sensor;  check  output  with  scope  and 
counter. 

4 

Check  level  shifter. 

5 

Check  period  multiplier. 

6 

Trim  reference  oscillator  to  ±0.  1 Hertz. 

7 

Check  scale  selector.  See  Figure  9. 

8 

Check  product  modulator.  See  Figure  11. 

9 

Check  subtractor.  .See  Figure  11. 

7.2  liiti-rfai'e  Tests 

Tests  described  in  this  section  were  accomplished  by  replacing  the  sensor 
output  with  a function  generator  monitored  by  a frequency  counter.  Application  of 
a very  gentle  touch  on  the  generator  vernier  control  and  a good  deal  of  patience, 
enabled  the  periods  to  be  set  to  five  place  accuracy. 

The  initial  data  were  gathered  based  on  the  periods  in  Table  2.  The  results 
were  repeatable  within  a bit,  but  above  12.  2 km  (40  k ft)  the  error  grew  to  -2  bits, 
above  T.T.  5 km  (110  k ft)  to  -10  bits,  and  above  .Tfi.  6 km  (120  k ft)  to  -20  bits. 
Analyzing  the  test  setup  led  to  two  conclusions:  (1)  at  the  highest  altitudes  the 
error  was  due  to  the  accuracy  of  the  input  period  which,  in  turn,  was  a function  of 
Table  2 and  the  function  generator;  and  (2)  in  the  middle  range  12.2  km  to  18.3  km 
(40  k ft  to  60  k ft),  the  error  was  due  to  the  linear  interpolation  error  when  Table  2 
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was  calculated  from  tables  in  the  U.  S.  Standard  Atmosphere  (reference  2)  and  the 
sensor  instruction  manual  (reference  1). 

The  function  generator  was  the  easier  of  the  two  problems  to  solve.  It  was 
replaced  with  a more  stable  unit  which  could  hold  better  than  six  places  accurately. 
Producing  more  accurate  periods  for  test  altitudes  took  a good  deal  of  time. 

Six  place  accuracy  was  the  minimum  acceptable.  The  U.  S.  Standard  Atmos- 
phere is  tabulated  to  five  places;  the  sensor  manual  has  data  to  six  places.  How- 
ever, since  neither  function  is  linear,  interpolation  would  introduce  an  error  due 
to  the  curvature  of  the  function.  Maximum  accuracy  could  be  obtained  only  by 
calculating  the  exact  period.  Calculating  atmospheric  pressure  at  some  altitude 
is  straightforward;  calculating  sensor  period  from  that  pressure  is  more  involved. 

The  sensor  function  is  defined  in  section  2.  To  proceed  from  the  dependent 
variable  (pressure)  to  the  independent  variable  (period),  a point  search  by  interval 
halving  algorithm  was  devised.  This  algorithm  produced  a period  which  corre- 
sponded to  a pressure  accurate  to  within  one-tenth  the  absolute  sensor  error.  The 
algorithm  flowchart  is  in  Appendix  E. 

The  standard  atmosphere  formula  and  the  sensor  formula  with  the  point  search 
algorithm  were  coded  into  a Hewlett-Packard  9810  calculator,  and  a new  period  vs 
altitude  table  was  produced.  Comparing  the  periods  in  Table  2 and  Table  13,  high- 
lights the  increased  accuracy  of  the  latter  data.  Table  13  holds  the  data  generated 
when  the  interface  was  then  retested  using  the  more  accurate  periods.  If  the 
empirical  data  are  compared  to  the  computer  printout  in  Appendix  C,  it  will  be 
noted  that  every  point  is  within  less  than  one  bit  of  the  control  altitude.  This  test 
effectively  closes  the  loop  on  the  interface  design. 

7.3  \ltimeler  tlhumber  Tests 

The  final  test  includes  the  sensor  and  interface  in  an  "altitude"  check  via 
vacuum  chamber.  This  test  was  accomplished  using  three  different  instruments 
to  monitor  the  chamber  pressure  altitude.  A Wallace  and  Tiernan  (W  & T)  Model 
FA129  gauge  was  used  from  1.  52  to  12.  2 km  (5  to  40  k ft);  a MKS  baratron  pres- 
sure meter  type  77,  from  12.  2 to  25.  9 km  (40  to  85  k ft);  and  a W & T Model 
FAIRO  gauge,  from  25.  9 to  44.  2 km  (85  to  145  k ft).  In  Table  14,  Altimeter  Test 
Data,  the  points  where  the  monitor  instrument  changed  are  listed  twice,  one 
datum  point  on  each  instrument.  The  reported  altitude  is  the  mean  of  data  gathered 
by  three  technicians.  In  all  cases,  the  data  are  within  the  absolute  error  of  the 
altimeter/chamber  system. 

Conspicuous  by  its  absence  is  a datum  point  at  0 km.  Recall  Table  7.  the 
sensor  accuracy  is  1.  5 m (5  ft)  at  that  altitude.  Also,  the  altimeter  resolution  is 
18  m (58  ft)  at  that  altitude.  Therefore,  at  sea  level  the  altimeter  readily  responds 
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to  changes  in  barometric  pressure.  An  ambient  check  of  its  accuracy  was  accom- 
plished using  a quartz  manometer.  In  every  case,  the  closure  was  less  than  the 
units  resolution. 


Table  13.  Interface  Test  Data 


Altitude 
k ft 

Period 
^ sec 

Empirical  data 
k ft 

0 

196.  0438 

0.02 

5 

201.3297 

4.  98 

10 

206. 3008 

10.  01 

15 

210.  9188 

14.  95 

20 

215.  1523 

19.  97 

25 

218.  9820 

24.  94 

30 

222.3988 

29.  96 

35 

225.4039 

35.  00 

40 

227.9750 

39.  99 

45 

230.0862 

45.02 

50 

231. 8056 

49.97 

55 

233. 1973 

54.  96 

60 

234.3181 

59.97 

65 

235. 2167 

64.  96 

70 

235.  9329 

p9.  98 

75 

236. 5000 

75.  00 

80 

236. 9489 

73.96 

85 

237. 3039 

85.  01 

90 

237. 5847 

90.02 

95 

237.  8069 

94.  99 

100 

237. 9829 

100.05 

105 

238.  1223 

105.  03 

TOO 

238. 2322 

110,01 

115 

238.3187 

115.04 

120 

238. 3871 

120.05 

125 

238. 4412 

125.  02 

130 

238.  4843 

130.00 

135 

238. 5188 

135.  04 

140 

238.  5464 

140.00 

145 

238. 5686 

145.  06 
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Table  14.  Altimeter  Chamber  Test  Data 


= 

Control  Altitude 

Reported 
Altitude 
k ft 

Standard 
Deviation 
k ft 

km 

k ft 

1.  52 

5 

5.02 

3.  05 

10 

9.  97 

4.  57 

15 

15.  04 

0.04 

f).  10 

20 

20.09 

0.  05 

7.  62 

25 

25.05 

0.05 

9.  14 

30 

.29.  96 

0.05 

10.  7 

35 

34.  91 

0.  03 

12.  2 

40 

39.88 

0.05 

12.  2 

40 

40.  05 

0.00 

13.  7 

45 

45.  07 

0.  05 

15.  2 

50 

50.04 

0.  03 

16.  8 

55 

55.  04 

0.  04 

18.  3 

60 

59.  97 

0.00 

19.  8 

65 

64.  89 

0.03 

21.3 

70 

69.  81 

0.03 

22.  9 

75 

74.  69 

0.00 

24.4 

80 

79.  60 

0.05 

25.  9 

85 

84.42 

0.03 

25.  9 

85 

84.  95 

0.  06 

27.  4 

90 

89.  90 

0.09 

29.  0 

95 

94.  77 

0.06 

30.  5 

100 

99.  76 

0.  13 

32.  0 

105 

104. 67 

0.05 

33.  5 

no 

109.  59 

0.06 

35.  1 

115 

114.  60 

0.  13 

36.  6 

120 

1 19.  50 

0.  28 

38.  1 

125 

124.  72 

0.  21 

39.  6 

130 

129.  28 

0.31 

41.  1 

135 

134.  79 

0.  79 

42.  7 

140 

139.  69 

1.  13 

44.  2 

145 

145.  60 

0.  59 
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Tables  13  and  14  are  the  best  conclusion  this  report  could  have;  the  instrument 
worked  as  designed.  Figure  15  depicts  the  functioning  of  the  altimeter  and  compu- 
ter program  plus  the  tests  documented  above.  At  present  only  one  of  these  instru- 
ments has  been  fabricated,  but  the  second  will  be  schematically  identical  to  the 
first.  The  firmware  (ROM's)  and  software  (altitude  dictionary)  will  change  in  re- 
sponse to  the  coefficients  of  the  sensor  equations  as  set  forth  in  the  algorithm. 


CONTROL  altitude  CONTROL 

ALTITUDE  DICTIONARY  ALTITUDE 


CHECK:  TABLE  13  CHECK:  TABLE  14 


Figure  15.  Altimeter  System  and  Tests 
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Appendix  A 

As  Built  Schemotics 


Figure  Al.  Altimeter  Block  Diagram 
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-15  VDC 
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73 


AMPHENOL 

223-1609 


RED 

GRN 

WHT 

BLK 


♦ HAMILTON  STANDARD 
KEY  = 

0 , CONNECTION  NUMBER  (CN) 
a , CN  QUANTITY 
- 2 EA 
y - 3EA 
8 -4EA 
w'-30EA 


@ UNUSED 

Q BOARD  CONNECTOR  PIN 


CHASSIS 


Figure  A2.  Sensor  and  Case  Connectors 


4r> 


BOARD  3 


Figure  A7a.  Subtracter  Control  and  Counter 
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SCALE:  1357  0246 

Figure  A7b.  Subtracter  ROM 
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TYP 


C56 
24  pf 


C57 


C58  C59 


C60  C6I 


C62  C63 


ROW 


X 


T 


0 

r 

) 13 

LXJ 

0 
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0 
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COL 
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X 


12 


X 


10 
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R23 

56000 
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W (M  M (M. 

XTXTxlxlxIxTjrlr' 
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CRI2 
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12 
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* NOTE:  0-3  A30:  CD40I6AE 

4-1  A3l:  CD40I6AE 

PIN  7 GNO 


14  V 


DO 


BOARD  3 


Figure  A7c.  Subtracter  Switch 
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SENSOR  INTERFACE 


ENCODER  DATA  CHANNEL 


Figure  A8.  Interiace/Encoder  Interconnection 


53 


Computer  Program  for  Dictionary 


I 

I^EDINO  PAai  BUNK-MOT  Wl 


TM6  INPUT  DATA  ARC  FOR  HAMILTON  STANDARD  DIGITAL  PR 


A*4  -tSA.  T ( INPUT.  OUTPUT) 


OMFNSIOM  Li  TTRIfi).  ISai  0)  « A15)  ,HHI9),HF(S)  fSI2) 
OAIA  LETT ^/IHS.IMU.IMR.IHH.IMD.IHK. IHG • IHO/ 


REAa  121,>>l.AiB 
■ FJRMAr<lAi0.5F10.5,2Fl0.5) 


eFAO  120, IS3 

120  F0-<HATCdI10l 

PRINT  TOO 
F0PN4T  I H 


irSSURE-) 
PRINT  T01,3'< 


TOl  FORHATIl^O.iK, 

istants-m 

P<INT  TO? 

FORMATURJ 


PRINT  T3r 

T30  Fj^NATTHO.eex  ,-Y=-B4  0)  ♦B(  1)  •PERIOD-) 


PRINT  Til 

FORMAT  I HO  ,bX,-AR£“» 


PRINT  T3?, 1) 

T3?  F JRRAT  CIMJ,  31X ,-0(Q) a-.Fld.S  ) 

PRINT  T13,3I2) 

Til  FORHATIHJ.JIK  , - E 1 1 1 s-,  F I 0 . 5 ) 


TPANSOUCCH  HOOtL  10 


FOR  TMC  SENSOR  CALIBRATION  EQUATION-) 


RINT  TOl.Ad) 

TJI  F ORMAT  ( Hj  ,31X  ,-As-,  FIO  .5) 


PRIST  Tai»,A<2l 
Tj4  FORMAT  ( HO, 31X  ,“B 


PRINT  Ti}p,A<3) 

TJ5  FJRMAT I IMO, JIX - - 


PRINT  TOf,A(A) 
Tgf»  FORMAT  T 1M0,31X  , -I 


PRINT  TOT. ATS) 
TOT  FORMATIIMO.SIX. 


PRINT  Til 

FORMAT(l'«Q.bX.-THE  SUBTRAHENOS  FOR  THE  EIGHT  SCALES  ARE“) 


PRINT  T1?,I>B<1) 

Ti2  F0RMAT<HQ.2bX. -scale  Ql-.IlO) 


PRINT  T13, 136(2) 

T13  FORMAT ( IMP. ?6X. -SCALE  1I-,I10) 


■#  U' TK4-l<nfrT*LlM. 


PRINT  T16,1S3<S) 

F 0 R M A T ( 1 H 0 , 2 8 X , -SCALE4J  - , 1 1 0 ) 


I 


pf^o&R4ii  hSALt  7k/Th  oAt*!  TTN  V.MWI  12/13/77  12.37A7 

PRINT  Tl?,I>8(a> 

jo 

719 

FORMAT(l’43«26Xi**SC*LC  7l**»lt0) 

PRINT  720 

720 

FORHAT<1HO,6X,**THE  CQUATION  AND  CONSTANTS  lMIC»4  DEFIME  THE  U.S. 

IASOARD-) 

721 

PRINT  721 

FORMAT(lHO«3Xt**ATNOSPHeRC,1976  TO  47  KN  ARC  EHSEOOEO  IN  THE  PROCRA 

66 

IM.'M 

I’GaO 

DO  iOQ  is:sl«b 

I J<=0 

ITC=IS3(I53> 

NaIiiC*3 

UD  600 

I-»&=IP3tl 

101 

PRINT  101, 3N, IPG 

F3RNAT(1H1,6X,**G£0P0TENTIAC  ALTITOOE  VERSUS  COD€*  ,OX,**SN*,  1 AlO  t OX, 

76 

l”oAGE”, 13, 20X, “AVERAGE  RESOLUTlblTI 

PMNT  1)2 

102 

FDRHAT  1 HO,30X,-ALTITUOEI  K1LOHETERS",50X,^HE7W5"1 

PRINT  lOF 

90 

105 
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PRINT  IJ3,L-TTR 

103 
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nj  400  J7sl,2 

I U»I0<»1 

P>INT  13«. 

J6 

1U4 

FORNATCH  » 

PRINT  l(jc 

1 J = 

FORNATCIH  ) 

33  304  J3*1,B 

0)  200  J4:  l,c 

P3sITC72.**^ 

y = pj*e(2)  -61 1> 
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2U1 
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35 
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204 
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Appendix  D 

A Method  of  Subtraction 

I.  TIIK  MKTHOI) 

The  circuit  in  Figure  Dl(a)  was  presented  to  the  author  by  Mr.  Hans  Laping 
■AFGL/LCC  as  a method  of  determining  the  difference  between  an  unknown  frequency, 
fu  (U  events/second),  and  some  fixed  standard,  f^  (N  events/second),  when  the 
unknown  is  larger  than  the  standard.  A 1 sec  sample  of  the  signal  of  unknown  fre- 
quency is  applied  to  a counter  which  addresses  an  AND  gate  (Gl).  This  gate  detects 
a count  equal  to  N events  and  sets  a latch  to  "remember"  the  event.  Gate  G2  does 
the  actual  frequency  comparison  by  passing  any  input  signal  which  occurs  after  the 
Nth  event  is  detected.  Two  types  of  outputs  are  available;  a pulse  of  variable  length 
if  switch  .S  is  in  the  "a"  position,  or  a number  of  pulses  (F-N)  if  S is  in  the  "b" 
position.  In  our  case,  the  latter  method  is  used  and  the  pulses  are  counted  as  men- 
tioned in  Section  1 of  the  basic  report. 
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that  has  the  advantage  of  programability.  This  lends  the  instrument  flexibility  to 
accept  various  (serial  number)  sensors  without  changing  wiring.  Each  sensor  has 
its  own  companion  ROM  which  is  placed  into  the  circuit  with  the  sensor. 


Appendix  E 

Pngrommoblt  Calculator  Flow  Chort 


Since  there  is  a large  number  of  programmable  calculators  available,  only  the 
Flow  Chart  of  the  fKjint  search  program  will  be  given  (see  Figure  El).  It  will  work 
for  any  monotonic  curve  up  to  the  fourth  order  and  can  handle  a dummy  variable  z 
where 

z = B + B,  X , 
o 1 ’ 


and 


y = + Aj  z + A2  + Ag  z^  + z^  . 

If  the  dummy  variable  z is  not  used,  then  z = x,  B^  = 0 and  Bj  = 1. 

Figure  E2  is  an  annotated  curve  for  the  application  under  discussion.  The 
subscripts  upper,  u,  and  lower,  1,  are  defined  in  terms  of  the  dependent  variable 
pressure,  P,  and  then  applied  to  the  corresponding  ordinate  period  p.  In  the  pro- 
gram, the  pressure  of  interest,  P^,  is  calculated  as  a function  of  the  standard 
atmosphere.  The  two  bounding  points  on  the  curve  are  chosen  from  the  sensor 
manual  and  entered  into  the  upper  and  lower  registers  for  each  variable.  The  pro- 
gram then  halves  the  interval  on  the  period  axis  until  the  P^  computed  from  a is 
within  the  limit:  e of  P^. 
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Figiare  Ela.  Flow  Chart 


Figure  Elb.  Flow  Chart 
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Figure  Elc.  Flow  Chart 
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